A Steady State Thermodynamic Entanglement Witness 
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We present a method for entanglement detection in steady state systems using a thermodynamic 
witness. To ihustrate, we consider an example and find that for an XX spin chain, the presence 
of an energy current increases the region of entanglement detected by the steady state witness. 
Further, we find that entanglement exists even at a high steady state temperature. We also discuss 
concurrence in the steady state system and find that the amount of entanglement can be increased 
by the current. 
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Entanglement is one of the most fascinating aspects 
of quantum physics and is an important resource in the 
field of quantum information. While a separable state 
can be written as a convex sum of pure product states, 
(J = ^^Picrj (8) cr| (g) • • • (g) (jj^, an entangled state can- 
not. Many-body entanglement [1] is not easy to quantify; 
while thermal pairwise entanglement can be measured us- 
ing concurrence, measuring multipartite entanglement is 
a difficult task. Entanglement witnesses are a useful al- 
ternative to these measures. An entanglement witness is 
an expectation value of an operator for which a bound 
can be found for any separable state. 

The subject of quantum thermodynamics (H has al- 
ready contributed to the success of entanglement wit- 
nesses [l, For example, magnetic susceptibility can 
witness entanglement [5], as can an average over the 
nonequilibrium work done during a process [6]. At 
present, thermodynamic witnesses are used to detect en- 
tanglement in equilibrium systems. Our aim is to ex- 
tend the success of equilibrium entanglement witnesses 
to nonequilibrium steady state systems. 

In this letter, we demonstrate how a steady state en- 
tanglement witness can be found. Specifically, we find a 
steady state thermodynamic quantity, a current, which 
can detect entanglement in steady state systems. We 
consider an XX spin chain in the thermodynamic limit 
to exemplify this witness and find that the introduction 
of an energy current increases the region of entanglement 
detected. In order to gain insight into whether the en- 
tanglement itself increases with the driving field, we also 
calculate the chain's correlation functions and use them 
to determine the nearest and next nearest neighbour con- 
currence. We discover that the concurrence can be in- 
creased by the presence of the energy current. 

Thermodynamics has been used in many areas of 
physics. It is limited to equilibrium systems, although 
some thermodynamic concepts such as work done apply 
during dynamical processes. No equivalent formulation 
for nonequilibrium thermodynamics exists at present. 
However, attempts have been made to provide a frame- 
work for steady state thermodynamics [3, [H since these 



systems are time independent and simpler to investigate 
than the general out of equilibrium case. There are also 
methods by which steady state systems can be modeled. 

In particular, the Hamiltonian of an equilibrium sys- 
tem can be modified to allow a steady state current pass- 
ing through the system to be described (9l-[ll|. In this 
way, nonequilibrium steady state systems can be mapped 
onto equilibrium systems. For an energy current, this 
is achieved by defining the operator of the current us- 
ing continuity equations, dh{x,t)/dt + dj{x^t)/dx = 
where h{x^t) is the energy density operator and j{x^t) 
is the heat flux operator. For a chain with N sites, the 
discrete form of this is dhi/dt = ji-i — ji^ and the left 
hand side can be written dhi/dt = i[H^hi\. For one di- 
mensional systems = + ^ ^^))-> ^^^^ 
been shown to be 



) no 



V{1 + 



lJ + 2)] + [V{l-lJ),V{l,l + l)]). (1) 



Here, is the non-interacting part of the Hamiltonian, 
and V{l^l + 1) is the interaction. The energy current 
operator is = "^pi and can be incorporated into 
the Hamiltonian formalism using a Lagrange multiplier, 
A, so that we have H = Hq — XJ^ . In the case that 
[J^^Hq] = 0, and Hq can be diagonalised analytically, 
the total Hamiltonian H can also be diagonalised using 
the same method. 

We note that the origin of the steady state current 
is not important when using this method. That is, the 
current is given by Eq. [T] whether it is induced by an 
external driving field or by a reservoir [l3[. For example, 
the steady state could be achieved by holding opposite 
ends of an open spin chain at different, constant, tem- 
peratures, but the formalism above is not limited to this 
scenario. 

Using Eq. [H the ground state correlation functions of 
an Ising spin chain have been investigated 13, [HI . Also 
at zero temperature, a steady state XX spin chain with 
an energy current increases the entanglement present [HI 
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compared to the equilibrium case, and quantum state 
transfer can be improved [l^- In this letter, we take 
this formulation of steady state systems into the thermal 
regime using the method outlined below. 

It is possible to construct a steady state thermal den- 
sity matrix 18, 19| which is similar in form to the equilib- 



rium case: p = e-/5^o-7^ jz^ where Z = tr(e"^^o"^'^ ) 
is the steady state partition function, is the steady 
state energy current operator discussed above, and /3 and 
7 are Lagrange multipliers. f3 = 1/T can be thought of 
as a generalised inverse temperature which is valid in 
steady state systems, while 7 = —X/3 is the driving term 
of the energy current and has no analogue in equilibrium 
thermodynamics. 

This definition of the steady state density matrix al- 
lows a nonequilibrium steady state entanglement witness 
of the form 



To aid in our later calculations for the witness, we 
introduce an extra term, 6, into the Hamiltonian such 
that Eq. [3]becomes Hxx = (^f ^f+i + ^f+i) ~ 

bB'^iCrf. The total Hamiltonian, can again 

be diagonalised using the transformations described 
above. The partition function is thus In Z = 
£ J^"" dq\n[2cosh{(3^{q))] where ^{q) = Bb- Jcosq^ 
2^JBsmq — 2^ sing cos Once the calculations are 
performed, we set 6 = 1 in order to recover the original 
system. 

We now introduce two entanglement witnesses. The 
first is similar to existing witnesses, [H, ^nd is found 
by calculating the expectation value of the Hamiltonian 
and rearranging the resulting equation. The second is 
found using Eq. [2] and is a genuine steady state witness. 
We consider both to demonstrate why the steady state 
witness is important. 

The first witness is 



(77 



(2) 



where ( and 77 are constants and Q = (J^) = tr( J^p) 
is the expectation value of the energy current, to be cal- 
culated. Thus we can detect entanglement in the steady 
state system using the energy current, a steady state 
quantity. We note that other currents, such as a mag- 
netisation current, can be calculated in a similar way to 
that described above. Therefore other currents can be 
used to witness entanglement similarly to Q in Eq. O 

To illustrate how the witness works, we consider an 
XX spin chain in the thermodynamic limit. 
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where J is the coupling strength between nearest 
neighbour sites and B is an external magnetic field. Us- 
ing Eq. [TJ the energy current for Hxx is 
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(4) 



We consider J and B to be positive throughout the 
letter. The total Hamiltonian, H = Hxx + -^J^ can be 
diagonalised [20| using a Jordan- Wigner transformation, 

— nm=i ^ {^1 ~ ^^^) /2 and a Fourier transform, 
-^Y^k^ke^^ to give H Yki^^kdldk - IB) 
where A{q) = (5 — J cos (2^ J sing + 1) as 

N ^ 00. Although we consider the thermodynamic 
limit, our method for calculating a steady state witness 
can be applied to a system of any size. 



Wi = 2 



U^BM-^Q 



JN 



(5) 



where U = {H) 
nal energy, U = - 



is the steady state equivalent of inter- 
^ In Z — ^ ^ In Z. The second term. 



5 ^ In Z is equivalent to and therefore cancels with 



term in Eq. [5l We also have — ^ In Z = 

^/o 6i(7(Jcosg-5)tanh(/3A((7)). We use the conven- 
tion that the magnetisation is given by M = Yi{cf^) = 
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InZ rather than by -^^InZ, though each 

definition gives the same witness on appropriate consid- 
eration of the remainder of the expectation value of the 
Hamiltonian. Thus M/N = J^"" ^ tanh [f3A{q)]. 

The system is entangled when Wi > 1. This bound 
is calculated using the expectation value, Kcrfaf^^i^ -h 
c^fcrf+i)|, which is equal to Eq. [5l The bound is found 
for pure product states [l, Q using the Cauchy-Schwarz 
inequality, i^iXiyi)'^ < Yixf^^yf and the definition 
of a single site density matrix which leads to the inequal- 
ity (erf )^ + (^z^)^ + (erf )^ < 1. Since the set of separable 
states is convex, this bound is also true for all separable 
states. 

Fig. [1] shows the region of entanglement detected by 
witness Wi. It shows that increasing the driving quan- 
tity 7 increases the region of entanglement detected by 
the witness. This is an interesting result since it indicates 
that introducing an energy current increases the entan- 
gled region in this thermal system. That is, 7 increases 
which values of 5, T etc. the system is entangled for. 
In addition, we see that entanglement can exist even at 
high steady state temperatures with a large enough 7. 
However, witness Wi is not easily experimentally mea- 
sured since U and M are now steady state thermody- 
namic quantities and thus no longer correspond to clear 
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FIG. 1: The surface in the plot corresponds to Wi = 1. The 
witness detects entanglement within the convex regions where 
Wi > 1. We consider J — 1 and T is the generalised steady 
state temperature. 



FIG. 2: The surface in the plot corresponds to Wss = 1. 
The witness detects entanglement within the convex regions 
where Wss > 1- We consider J = 1 and T is the generalised 
steady state temperature. 



measurable quantities. Therefore, we use Eq. [2] to calcu- 
late a genuine steady state entanglement witness, 



JN{2B + J) • 



(6) 



Here we use the average energy current, Q, itself to 
detect entanglement. Since Q is the expectation value 
of Eq. O we can again use the Cauchy-Schwarz in- 
equality and the definition of the density matrix to 
find a bound for separable states. This leaves \Q\ < 
BJN-\-^ I (erf ) I for pure product states. For any pure 
product state, |(crf)| < 1, hence we can rearrange the 
inequality and find that entanglement can be detected 
when Wss > 1- The expectation value of the energy 
current is Q = ^ J^^ dq{ J cos q — B) smqtdinh.{f3A{q)). 
The steady state witness is shown in Fig. [2l Again in- 
creasing 7 increases the region of entanglement detected, 
and entanglement is detected even at high steady state 
temperatures. For this witness, no region is detected at 
low values of B and T, but an extra region exists when 
B is close to zero and higher T. An interesting way to 
view this witness is to consider Wss itself. It shows a 
large, absolute, expectation value of the energy current 
will detect entanglement (after scaling by the constant 
r] = 2/[JN{2B^J)]). 

How the expectation value of the energy current would 
be measured experimentally is dependent on the cause 
of the steady state. For example, for a spin chain in a 
constant temperature gradient, the heat conduction and 
hence Q can be experimentally determined [21 1. 

By considering the concurrence, we can confirm that 
it is the amount of entanglement which increases with 7. 



Concurrence quantifies entanglement between two mixed 
qubits, and is given by C{p) = max{0, Ai — A2 — A3 — A4} 
where Ai > A2 > A3 > A4. The A^s are the square roots 
of the eigenvalues of the matrix pp where p = {a^ 
cry)p*{cFy (g)cr^). Since [H.^icrf] = 0, the concurrence 
is given by C{pi^i^r) = 2max{|z| - y^, 0} ^ where 



I "^l^R 



))/4, vy : 



((1 + {crfcrf^ji))'^ - 4(crf)2)/16, and we have used that 
(<jf ) = {crf_^i) due to the chain's translational invar iance. 

We calculate each of the correlation func- 
tions following the method in [23], and calculate 
the nearest neighbour, R = 1, and next near- 
est neighbour, R = 2, concurrence. Generally, 

^r^i^) = -^e^^^'iAiU^ZliAi^rnBi^M^n) + 
^e'^^'iBi Ul~JMi^mBi^m)Bi^R). {cTfcTf^n) 
(AiBiAi^rBi^r) and (af ) = (AiBi) where Ai = {aj^ai) 
and Bi = {aj — ai). Using Wick's theorem, we can 
rewrite each of these in terms of two point correlation 
functions. For example, the zz correlation function 
is {AiBiAij^rBij^r) = - G| + S'l where we have 
defined Gr = {AiBij^r) with Gr = -{BiAij^r) and 
Sr = {BiBi+r) = -{AiAi^r). 

Using the same method of diagonalisation as for the 
Hamiltonian, we can write (af) = ;^X1/.(1 — 2d\dk)- 
In addition, Gr = ^EfcCos(^)(l - 244) and 
Sr= EfeSin(^)(l - 2dldk). Therefore, the ther- 
modynamic expressions for Gr and Sr are found directly 
from the magnetisation, M/N = (af ) calculated previ- 
ously: 
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FIG. 3: Nearest and next nearest neighbour concurrence re- 
spectively with a) and d) when 7 = 0, b) and e) when 7 = 1, 
and c) and f) when 7 = 2. We have set J = 1 throughout. 



27r 



dq 



Gr = I ^cos(gi?)tanh(/3A) (7) 

Sr = i f — sin(gi?)tanh(/3A). (8) 
Jo 27r 

We plot the concurrence against the magnetic field and 
steady state temperature in Fig. [3l We find for nearest 
neighbour concurrence, Figs. [3] a), b) and c), that the 
amount of entanglement decreases at low T and B with 
increasing 7, and also that the region of entanglement 
decreases. The amount and region of entanglement in- 
creases with 7 when T and B are higher. The plot co- 
incides with both witnesses which detect entanglement 
at high B and T. Fig. [T] shows that the witness also 
detects a region of entanglement at low B and T that 
decreases with increasing 7 which is reflected in Fig. [3l 
Fig. [2] however, also detects a region of entanglement 
at low magnetic field, and higher T and 7, which indi- 
cates this witness detects entanglement which cannot be 
categorised as nearest neighbour. Therefore we consider 
next nearest neighbour concurrence plotted in Fig. [3]d), 
e) and f). Both the amount and region of entanglement 



are reduced compared to nearest neighbour concurrence. 
The peak at low temperature remains constant for each 
value of 7, while as 7 is increased, new regions of entan- 
glement appear at low temperature and higher magnetic 
field. Thus for next nearest neighbour entanglement, 
the introduction of an energy current again increases the 
amount of entanglement. However, this does not explain 
the region of entanglement at low B and higher T and 7 
detected in Fig. O It would be interesting to determine 
the type of entanglement here. 

We have introduced a nonequilibrium steady state en- 
tanglement witness and using an XX spin chain, have 
demonstrated how such a witness works. We have shown 
that introducing an energy current increases the entan- 
glement detected in this system at high steady state tem- 
perature, and increases entanglement itself on consider- 
ation of the concurrence. 
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